Mon. Not. R. Astron. Soc. 000, 1-10 (2009) Printed 26 February 2010 



(MN MpjX style file v2.2) 



The rotational excitation of HCN and HNC by He: New 
insights on the HCN/HNC abundance ratio in molecular 
clouds 

E. Sarrasin 1 , D. Ben Abdallah 2 , M. Wernli 3 , A. Faure 3 , J. Cernicharo 4 and F. Lique 1 ' 5 ^ 

1 LOMC - FRE 3102, CNRS- Universite du Havre, 25 rue Philippe Lebon, BP 540, 76058, he Havre, France 

2 Laboratoire de Spectroscopic atomique, moleculaire et applications, Faculte des Sciences Universite Tunis el Manar, Tunis 1060, Tunisie 
3 Laboratoire d'Astrophysique de Grenoble, Universite Joseph Fourier, CNRS UMR 5571, BP 53, 3804-1 Grenoble Cedex 09, France 
4 Department of Astrophysics. CAB, INTA-CSIC. Crta Torrejon a Ajalvir km 4- 28850 Torrejon de Ardoz. Spain 
5 LERMA and UMR 8112 of CNRS, Observatoire de Paris-Meudon, 92195 Meudon Cedex, France 

Accepted 2009 month XX. Received 2009 month 06; in original form 2009 month 06 

ABSTRACT 

Modeling of molecular emission from interstellar clouds requires the calculation of 
rates for excitation by collisions with the most abundant species. The present paper 
focuses on the calculation of rate coefficients for rotational excitation of the HCN and 
HNC molecules in their ground vibrational state in collision with He. The calculations 
are based on new two-dimensional potential energy surfaces obtained from highly 
correlated ab initio calculations. Calculations of pure rotational (de)excitation cross 
sections of HCN and HNC by He were performed using the essentially exact close- 
coupling method. Cross sections for transitions among the 8 first rotational levels of 
HCN and HNC were calculated for kinetic energies up to 1000 cm . These cross 
sections were used to determine collisional rate constants for temperatures ranging 
from 5 K to 100 K. A propensity for even Aj transitions is observed in the case of 
HCN-He collisions whereas a propensity for odd Aj transitions is observed in the 
case of HNC-He collisions. The consequences for astrophysical models are evaluated 
and it is shown that the use of HCN rate coefficients to interpret HNC observations 
can lead to significant inaccuracies in the determination of the HNC abundance, in 
particular in cold dark clouds for which the new HNC rates show that the j = 1 — 
line of this species will be more easily excited by collisions than HCN. An important 
result of the new HNC-He rates is that the HNC/HCN abundance ratio derived from 
observations in cold clouds has to be revised from >1 to ~1, in good agreement with 
detailed chemical models available in the literature. 
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1 INTRODUCTION 

The isomers hydrogen cyanide (HCN) and hydrogen iso- 
cyanide (HNC) are among the most abundant organic 
molecules in space, from dark cold clouds (e.g. Irvine & 
Schloerb 1984) to circumstellar envelopes (e.g. Cernicharo 
et al. 1996), cool carbon stars (Harris et al. 2003), comets 
(e.g. Lis et al. 1997) and active galaxies (Perez-Beaupuits 
et al. 2007). HCN and HNC also belong to the small 
class of molecules detected in high-redshift galaxies, along 
with CO, CN and HCO+ (Guelin et al. 2007). In addition 
to thermal emission from various rotational transitions at 
(sub)millimeter and far-infrared wavelengths, a few maser- 
ing lines have also been detected for both isomers (e.g. Lu- 
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cas & Cernicharo 1989; Aalto et al. 2009). Thus, in contrast 
to purely thermochemical considerations (HNC is less sta- 
ble than HCN by about 0.6 eV, Bowman et al. 1993), the 
abundance of HNC in space is large: the abundance ratio 
[HNC]/[HCN] is observed to vary from ~ 1/100, e.g. in high 
mass star forming regions (Schilke et al. 1992), up to 4.5 in 
dark cold clouds (Hirota et al. 1998). The large variation of 
this ratio remains a puzzle in astrochemistry, although the 
major source of both isomers is thought to be the dissocia- 
tive recombination of HCNH + with electrons (e.g. Amano 
et al. 2008, and references therein). Finally, we note that the 
HNC/HCN (rotational) line ratio is now employed to dis- 
tinguish photon-dominated regions (PDR) and X-ray dom- 
inated regions (XDR): PDR sources all have ratios lower 
than unity while XDR have ratios larger than 1 (Meijerink 
& Spaans 2005). 
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Due to the low density of the interstellar medium (~ 
10 5 cm -3 ), the rotational levels of molecules are generally 
not at local thermodynamic equilibrium (LTE) . Reliable de- 
termination of molecular abundances therefore relies on ac- 
curate molecular collisional rates. Hydrogen molecules are 
generally the most abundant colliding parners in interstel- 
lar space, although collisions with H, He and free electrons 
can also play important roles in energetic regions. Rate co- 
efficients for the rotational excitation of HCN by He atoms 
(employed as substitutes for H2) have been calculated by 
Green & Thaddeus (1974), Green (unpublished) 1 and Mon- 
teiro & Stutzki (1986) (see also references therein). All these 
studies were based on the same potential energy surface 
(PES) obtained using the uniform electron gas model (Green 
& Thaddeus 1974). For HNC, there is to the best of our 
knowledge no collisional data except those of Faure et al. 
(2007) for electron-impact excitation. In astrophysical ap- 
plications, HNC and HCN have been therefore assumed to 
present similar collisional rates (see, in particular, the dis- 
cussion in Guelin et al. 2007). Similar HCN and HNC ro- 
tational rates were also observed for electron-impact excita- 
tion which is, however, dominated by the dipole interaction 
(Faure et al. 2007). In contrast, inelastic collisions between 
neutral species are dominated by short-range interactions 
and larger differences between HCN and HNC are expected. 

In the present work we present new rotational rate 
coefficients for HCN and HNC based on highly accurate, 
HCN— He and HNC— He PES. The paper is organized as 
follows: Section 2 describes the PES used in this work. Sec- 
tion 3 then contains a rapid description of the scattering 
calculations. In Section 4 we present and discuss our results. 
Finally, we analyze in section 5 the effect of these new rate 
coefficients on the excitation of HCN and HNC by modeling 
this excitation through a large velocity gradient code. 



2 POTENTIAL ENERGY SURFACE 
2.1 HCN— He 

Several theoretical PES for the HCN-He system have been 
published during the last fifteen years: Atkins & Hutson 
(1996); Drucker et al. (1995) and Toczylowski et al. (2001), 
hereafter called TDC01. In the present study, we used the 
latter surface. As a reminder, it was computed ab initio us- 
ing the CCSD(T) method (Hampel et al. 1992), with a triple 
zeta (aug-cc-pVTZ) basis set and an additional [3s3p2d2flg] 
set of bond functions (full details on the ab initio proce- 
dure used can be found in TDC01). The HCN molecule 
was treated as a linear rigid rotor with intramolecular dis- 
tances fixed at the equilibrium values thc='2. 01350 bohr 
and rcjv=2. 17923 bohr (see discussion below in the case of 
HNC). The global minimum of this PES is in the linear He- 
HCN configuration at J?=7.97 bohr and has a well depth of 
29.90 cm~\ 

1 The earlier calculations of Green & Thaddeus (1974) (re- 
stricted to the lowest 8 levels of HCN and temperatures of 
5—100 K) were extended in 1993 to obtain rate constants 
among the lowest 30 rotational levels and for temperatures 
of 100—1200 K. These unpublished results are available at 
http: //data. giss .nasa.gov/mcrat es/#hcn 



A semi-empirical surface (Harada et al. 2002) is the lat- 
est PES published for HCN-He. It used the TDC01 surface 
as a starting point. The choice of the TDC01 surface seemed 
more adequate to us, as it is purely ab initio and thus con- 
sistent in the whole relevant interaction range, while the 
empirical PES used the measured spectrum of the bound 
van der Waals complex to refine the TDC01 surface in the 
region of the attractive well, thus possibly losing precision 
in other regions that could be as important for the present 
dynamical applications. A contour plot of the potential is 
shown in Fig. 1 of TDC01. 

The PES being strongly anisotropic due to the length of 
the HCN rod, we fitted the TDC01 surface with a Legendre 
polynomials basis 

v(R,e) = J2 v *(R)Px(™se) (i) 

A 

using exactly the same truncation technique as for the 
HCsN-He system in Wernli et al. (2007a, b). In this equa- 
tion, R is the distance from the HCN centre of mass to 
the He atom, 9 being the He-HCN angle as measured from 
the H side of HCN. The combination of 25 polynomials 
(Amaa; = 24) and a cubic spline interpolation for the ra- 
dial coefficients (V\(R)) was employed to obtain a precise 
fitting of the PES, the standard deviation between fitted 
and ab initio data being < 1 cm -1 for all potential values 
lower than 300 cm -1 . We also tested our fit against the ab 
initio values given in the lower part of table I of TDC01, 
and found a difference ^ 3 cm -1 for all 16 points in this 
table. 

As a first outcome of the fitting, the V\ with even A 
were found to be significantly lower than those of Green & 
Thaddeus (1974) (table 3 in their article), indicating a lower 
(even) anisotropy of the TDC01 surface. This is expected to 
have an effect in the propensity rules, as shall be further 
discussed in section 4. 

2.2 HNC-He 

To our knowledge, no PES exists for the HNC-He van der 
Waals system. The ground electronic state of the weakly 
bound HNC-He system is a singlet with A' refection sym- 
metry. Within the ground electronic state, the equilibrium 
geometry of the HNC molecule is linear. Therefore, HNC 
can be considered as a linear rigid rotor. The HNC-He 
"rigid rotor" PES is described by the two Jacobi coordi- 
nates R, the distance from the centre of mass of HNC to 
the He atom, and 8, the angle between R and the HNC 
bond axis f, with 9 — corresponding to colinear HeCNH. 
The HN and NC bond distances thn and rjvc were frozen at 
their experimental equilibrium value thn — 1.8813 bohr and 
r NC = 2.2103 bohr (Huber & Herzberg 1979). As demon- 
strated by Lique & Spielfiedel (2007) for the CS-He system, 
for non-hydride diatomic molecules two-dimensional PES 
calculated for a frozen bond distance or obtained from full 
3D PES by averaging over the intermolecular ground state 
vibrational wavefunction are very similar. Consequently, in 
the present case we anticipate that restricting thn and tnc 
to their equilibrium value will introduce little error into the 
calculated inelastic rate coefficients. 

The potential energy surface (PES) was calculated in 
the supermolecular approach by means of a single- and 
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Figure 1. Contour plot of the 2D HNC-He PES as a function of 
R and 8. The energies are in cm -1 with the zero of energy taken 
to be HNC+He at infinite separation. 



double-excitation coupled-cluster method (Hampel et al. 
1992) (CCSD). In addition, we included perturbative con- 
tributions from connected triple excitations computed as 
defined by Watts et al. (1993) [CCSD(T)]. The calcula- 
tions were done with the MOLPRO 2006.1 package. For all 
three atoms we used the standard correlation-consistent po- 
larized valence-quadruple-zeta basis sets of Dunning (Dun- 
ning 1989) (cc-pVQZ) augmented with the diffuse functions 
of s, p, d, f and g symmetries by Kendall et al. (1992) 
(aug-cc-pVQZ). This basis set was further augmented by 
the [3s3p2ri2/l(7] bond functions optimized by Cybulski & 
Toczylowski (1999) and placed at mid-distance between the 
He atom and the HNC centre of mass. 

At all geometries the Boys & Bernardi (1970) counter- 
poise procedure is used to correct for basis set superposition 
error (BSSE). In this procedure the interaction energy is 
defined by 

V(R, 9) = Ehnc-hc(R, 9) - Ehnc(R, 9) - E Hc (R, 9) (2) 

where the energies of the HNC and He subsystems are deter- 
mined with the full (four atoms plus bond functions) basis 
set. 

Interaction energies at a total of 446 geometries were 
computed. The values of the radial scattering coordinate R 
ranged from 4 bohr to 25 bohr. The angular grid was uniform 
with a 15 degree spacing from to 180°. A contour plot of 
the potential is shown in Fig. 1. For this weakly-bound sys- 
tem the global minimum in the interaction energy was found 
to be -46.83 cm" 1 (7? = 7.30 bohr, 9 = 180°) corresponding 
to collinear He-HNC. We note the strong anisotropy of the 
PES. 

The calculated interaction energies were fitted by means 
of the procedure described by Werner et al. (1989) for the 
CN-He system. The largest deviations between the fit and 
the ab initio points occur primarily in the repulsive region of 
the PES. Over the entire grid, the mean relative difference 
between the analytic fit and the ab initio calculations is 
0.5%. The dependence of the potential energy surface on the 
He-HNC angle was fitted by the usual Legendre expansion, 
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Figure 2. Plot of first four radial Legendre expansion coefficients 
(A = 0...3) as a function of R. Solid lines denote HCN-He while 
dashed lines denote HNC-He. 



Eq. (1). From an ab initio grid containing 13 values of 9, we 
were able to include terms up to A max = 12. 

The dependence on R of the dominant coefficients for 
the HCN-He and HNC-He PESs is shown in Fig. 2. We ob- 
serve that, for the HCN-He system, the largest (in magni- 
tude) of the anisotropic terms (A > 0) corresponds to A = 2. 
This implies that, to a first approximation, the PES is sym- 
metric with respect to 9 — > n — 9. At the opposite, for the 
HNC-He system, the largest of the anisotropic terms corre- 
sponds to A = 1, 3 in the relevant range of R, that is between 
7 and 8 bohrs. This reflects a large (odd) anisotropy for the 
PES. Consequences for the rotational excitation will be dis- 
cussed in the next paragraph. 



3 DYNAMICAL CALCULATIONS 

The main focus of this paper is the use of the two fitted 
HCN-He and HNC-He PES to determine rotational exci- 
tation and de-excitation cross sections of HCN and HNC 
molecules by He atoms. In the following, the rotational quan- 
tum number of HCN and HNC will be denoted j. Despite the 
presence of a hyperfine structure in both molecule (Garvey 
& de Lucia 1974), we consider only the rotational structure 
of HCN and HNC. 

To determine the inelastic cross sections we used the 
exact close coupling (CC) approach of Arthurs & Dalgarno 
(1960). The integral cross sections are obtained by summing 
the partial cross sections over a sufficiently large number of 
values of the total angular momentum J until convergence is 
reached. The standard time-independent coupled scattering 
equations were solved using the MOLSCAT code (Hutson & 
Green 1994). Calculations were carried out at values of the 
total energy ranging from 3.2 to 1000 cm -1 . The integra- 
tion parameters were chosen to ensure convergence of the 
cross sections. We extended the rotational basis sufficiently 
to ensure convergence of the inelastic cross sections. At the 
largest total energy considered (1000 cm -1 ) the rotational 
basis was extended to j — 14 and 16, respectively, for the 
HCN-He and HNC-He calculations. The maximum value 
of the total angular momentum J used in the calculations 
was set large enough that the inelastic cross sections were 
converged to within 0.005 A 2 . 
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From the rotationally inelastic cross sections o" ■;_>-,' (E c ), 
one can obtain the corresponding thermal rate coefficients 
at temperature Tk by an average over the collision energy 
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where /i is the reduced mass and ks is Boltzmann's con- 
stant. To obtain precise values for the rate constants, the 
energy grid was chosen to be sufficiently fine to include the 
numerous scattering resonances which will be described be- 
low. 



4 RESULTS 

4.1 Cross sections 

Figure 3 illustrates the typical energy dependence of the 
collisional de-excitation cross sections obtained from the 
present CC calculations for a few selected rotational levels. 

The de-excitation cross sections are almost decreas- 
ing functions of the energy. For collision energies below 50 
cm -1 , many resonances are found. These are a consequence 
of the quasibound states arising from tunneling through the 
centrifugal energy barrier (shape resonances), or from the 
presence of an attractive potential well that allows the He 
atom to be temporarily trapped into the well and hence 
quasibound states to be formed (Feshbach resonances) be- 
fore the complex dissociates (Smith et al. 1979). Because 
of the averaging over collision energy, Eq. (3), these narrow 
resonances will have little, if any, effect on the relaxation 
rate coefficients. 
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4.2 Rate coefficients 

We obtained, by energy averaging, de-excitation rate coef- 
ficients for the first 8 (.7 =0 — 7) rotational levels of HCN 
and HNC, from the CC cross sections. The representative 
variation with temperature is illustrated in Fig 4. 

This complete set of (de)excitation rate coefficients is 
available on-line from the LAMDA 2 and BASECOL 3 web- 
site. Excitation rate coefficients can be easily obtained by 
detailed balance: 

= kj ,^(T K )^±^ex P [- ( Sj , - ej ) /kT K ] (4) 

where Ej and e 3 -/ are, respectively, the energies of the rota- 
tional levels j and j' . 

As they are the only complete set of rotational rates 
available in the literature, astronomers still generally use 
the HCN-He rates of Green & Thaddeus (1974) in order to 
interpret HCN and HNC observations. Table 1 compares on 
a small sample the temperature dependence of our HCN-He 
rates versus those of Green & Thaddeus (1974). We see that 
logically, the largest differences appear at low temperature. 

2 http://www.strw.leidenuniv.nl/ moldata/ 

3 http://basecol.obspm.fr/ 



Figure 3. Typical rotational de-excitation cross sections for the 
HCN [upper panel] and HNC [lower panel] molecules in collision 
with He as a function of the collision energy. 



These differences can be attributed to the use of a new PES 
for the scattering calculations. The maximum difference at 
100 K - the highest temperature for our data - is found on 
the smallest rate and is of a factor of 2. On the larger rates 
the two sets match nearly. 

4.3 Comparison between HCN and HNC rate 
coefficients 

Fig. 5 shows the HCN-He and HNC-He de-excitation rate 
coefficients from the j = 5 level at 10 K and 50 K. 

One can see that significant differences exist between 
the HCN-He and HNC-He rate coefficients: 

(i) On the overall, the HNC rate coefficients seems to be 
larger than the HCN rate coefficients. This can be explained 
by a larger well depth of the HNC-He PES than the one of 
the HCN-He PES (46.83 versus 29.90 cm" 1 ). 

(ii) The propensity rules seen in the two sets of rate 
coefficients are different. The HCN-He rate coefficients 
present a strong propensity in favour of transitions with 
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Figure 4. Typical rate coefficients for the HCN [upper panel] and 
HNC [lower panel] molecules in collision with He as a function of 
the temperature. 



even Aj whereas the HNC-He rate coefficients present a 
strong propensity in favor of transitions with odd Aj. 
These propensity are due to the shape of the PESs. Indeed, 
near-homonuclear symmetry of the potential energy surface, 
such as HCN-He PES, will favor transitions with even Aj 
whereas anisotropic potential energy surface, such as the 
HNC-He PES [see Fig. (2)], will favor transitions with odd 
Aj (McCurdy & Miller 1977). 

From this comparison, one can see that the use of HCN 
rate coefficients, in astrophysical applications, in order to 
interpret HNC observations may be dangerous since the use 
of "real HNC rate coefficients" will probably significantly 
modify the excitation of the HNC molecule. This will be 
discussed in the next part. 
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Table 1. HCN-He and HNC-He rate coefficients at 10 K and 
100 K in 10 -10 cm 3 molecule -1 s" 1 
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Figure 5. HCN-He and HNC-He de-excitation rate coefficients 
from the initial level j = 5 at 10 K [upper panel] and 50 K [lower 
panel] . 



5 ASTROPHYSICAL APPLICATIONS 
5.1 HCN 

We have estimated the line intensities of HCN using the 
present rate coefficients and compared the results with those 
obtained using the one of Green & Thaddeus (1974). Both 
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sets of rates have been corrected for the different reduced 
mass of the system HCN-H2, as molecular hydrogen is the 
main collision partner in molecular clouds : 

, , I MHCN-He \ 3 

«HCN-H 2 = KHCN-He (o) 

\ v A t HCN-H 2 / 

where /^hcn— He and aihcn-h 2 are respectivelly the reduced 
mass of the HCN-He and HCN-H2 collisional systems. 

We have used a Large Velocity Gradient (LVG) code 
to derive excitation temperatures, opacities, and brightness 
temperatures under different physical conditions. We have 
assumed a spherical cloud with an intrinsic linewidth of 
1 kms . We have selected three different temperatures: 
10, 30 and 50 K, and the density has been varied from 10 
to 10 9 cm -3 . The column densities of HCN had gone from 
10 12 (optically thin case) to 10 15 cm -2 (very optically thick 
case). The results using the new rate coefficients are shown 
in Fig. (6) (solid lines; dashed lines are the results for HNC 
- see next section -) which can be used to estimate physical 
conditions of the gas from the observed intensities of the 
,7 = 1 — 0, 2 — 1, and 3 — 2 lines in astronomical objects. 

In the optically thin case (iV(HCN) ^ 10 13 cm" 2 ) very 
high densities are needed to excite significantly the j = 1 — 
line of HCN while in the optically thick case, the HCN exci- 
tation is dominated by radiative effects and strong emission 
could be obtained even in the case of densities 10 4 cm" 3 . 
For these large opacities and narrow lines more sophisti- 
cated models taking into account the hyperfine structure of 
the molecule have to be performed to deal in detail with 
the radiative transfer of HCN (see Kwan & Scoville (1975); 
Guilloteau & Baudry (1981); Cernicharo et al. (1984); Cer- 
nicharo & Guelin (1987); Gonzalez- Alfonso & Cernicharo 
(1993)). 

In order to compare our results with those obtained us- 
ing the Green & Thaddeus (1974) rates, we show in Fig. (7) 
the ratio of the intensity for the j = 1 — 0, j = 2 — 1 and j = 
3 — 2 lines obtained using the new rate coefficients and those 
derived using the previous ones, R = Ts[New]/TB[01d]. 
Similarly to the calculations described above the number of 
rotational levels included in the calculations is 8 (j = — 7). 
The temperature varies between 10 and 100 K, the high- 
est temperature of the new and old HCN-He rate coeffi- 
cients. We have considered three different column densities 
for HCN, 10 12 , 10 13 , and 10 14 cm -2 , which correspond for 
the j = 1 — transition to the optically thin (r < 1), moder- 
ate optically thick (r ~ 1), and optically thick (r >1) cases 
respectively For the density we have selected also three val- 
ues, 10 4 , 10 5 and 10 6 cm , corresponding to low, moderate 
and high collisional excitation respectively. As one could ex- 
pect from a comparison of the rates given in Table 1, only 
small differences are found in the predicted intensity ratios. 
For the j = 1 — line we found similar intensities using 
both set of collisional rates at low temperatures (a few per- 
cent difference for T K < 20 K). The column density of HCN 
plays little effect in the low and moderate density cases. 
The largest differences for the j = 1 — line are found for 
n(H2) = 10 6 cm -3 . The situation is slightly different for the 
j = 2 — 1 and j = 3 — 2 lines. For low temperature the inten- 
sity ratios obtained with both set of rate coefficients show 
differences of 20 and 35% respectively for low and moderate 
densities (10 4 and 10 s cm" 3 ). 
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Figure 6. Solid lines: Brightness temperature for the j ' = 1 — 
0, j = 2-1, and j = 3-2 lines of HCN (indicated at the 
top of each column) derived using our new rate coefficients for 
temperatures of 10, 30 and 50 K (indicated in the left panels). 
The column density varies between 10 12 and 10 15 cm by steps 
of 10 1 / 2 - see text. Brightness temperatures are in K and densities 
in cm" 3 . Dashed lines: The same for HNC using the HNC-Hc 
rate coefficients obtained in this work. 



To further explore these effects we have considered 
the following case: T K = 50 K, n(H 2 ) = 10 5 cm" 3 and 
iV(HCN) = 10 13 cm" 2 . The old rate coefficients produce 
intensities for the j = 1 — 0, j = 2 — 1, and j = 3 — 2 
lines of 3.6, 2.2 and 0.8 K respectively. These intensities 
can be reproduced using the new rate coefficients for a den- 
sity of 1.3 10 5 cm" 3 . Hence, very small errors are intro- 
duced in the interpretation of the data using the Green & 
Thaddeus (1974) rate coefficients for warm molecular clouds. 
However, for Tk = 10 K the predicted intensities for these 
lines using the old rates is 1.56, 0.57 and 0.09 K respectively. 
With the new rates they are 1.55, 0.45 and 0.057 K respec- 
tively. The line intensity ratios j = 1 — 0/j = 2 — 1 and 
j — 1 — 0/j = 3 — 2 are different between both sets of rates 
while the j = 2 — 1/j = 3 — 2 intensity ratio is similar. Hence, 
the interpretation of astronomical observations of HCN us- 
ing the old HCN-He rates will be difficult in some cases as 
it will be hard to fit the observed line intensity ratios. Con- 
sequently, we recommend the use of the new HCN-He rate 
coefficients in interpreting observed line intensities. 



5.2 HNC 

Rate coefficients for HNC with Helium (or H2) have not 
been available until the present work. Interpretation of as- 
tronomical observations has been done using the HCN-He 
rate coefficients on the basis that the HNC-He and HCN-He 
PES could not be very different, that the dipole moment of 
both molecules differ by less than 3%, and that the energy 
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Figure 7. Ratio of the intensities of the j = 1 — 0, j = 2 — 1, and 
j = 3 — 2 lines of HCN (as indicated in the left panels) using our 
new collisional HCN-He rate coefficients over those obtained from 
the rates of Green & Thaddeus (1974). Three densities, 10 4 , 10 5 , 
and 10 6 cm -3 have been considered as indicated in the top panels. 
We have considered three cases for the opacity of the j = 1 — 
line : optically thin (TV(HCN) = 10 12 cm" 2 ), moderately thick 
(r ~ 1, AT(HCN) = 10 13 cm" 2 ), and optically thick (TV(HCN) = 
10 14 cm -2 ). The kinetic temperature varies between 10 and 100 K 
(see text). 



Figure 8. Ratio of the line intensities of the j = 1 — 0, j = 2 — 1, 
and .7=3 — 2 transitions of HNC (as indicated in the left panels) 
obtained using our new collisional rates HNC-He over those ob- 
tained when the HCN-He rates are adopted for HNC. Three den- 
sities, 10 4 , 10 5 , and 10 6 cm 3 have been considered as indicated 
in the top panels. We have considered three cases for the opacity 
of the j' = l — line at low kinetic temperature: optically thin 
(JV(HNC) = 10 12 cm" 2 ), moderately thick (r ~ 1, AT (HNC) = 
10 13 cm" 2 ), and optically thick (A 7 (HNC) = 10 14 cm" 2 ). The 
kinetic temperature varies between 10 and 100 K (see text). 



levels for both species are similar. However, the comparison 
of both set of rates done in previous sections indicate that 
some excitation effects could arise from the different propen- 
sity rules for both species and from the value of the feo->i 
rate coefficients (see Table 1), which differs by a factor ~ 5 
between both species. We have derived brightness tempera- 
tures for the HNC j = 1 — 0, 2 — 1, and 3 — 2 lines using the 
rates calculated in this work for the same set of models than 
for HCN. The results are shown Fig. 6 (dashed lines). The 
excitation temperature of the j = 1 — line of HNC is larger 
than that of HCN for the three temperatures considered in 
our calculations and for all column densities (abundances). 
The same applies to the j = 2 — 1 and j = 3 — 2 transitions 
for low temperature. However, the differences in predicted 
brightness temperatures for these two lines are less impor- 
tant for kinetic temperatures above 30 K. Hence, we have 
two different effects for HNC with respect to HCN: 

(i) the j = 1—0 line will be, for identical physical 
conditions and molecular abundances, considerably stronger 
for HNC than for HCN, i.e., it is more easily excited which 
is a direct result of the larger collisional rates for HNC. The 
same applies to the j = 2 — 1 and j = 3 — 2 lines but the 
effect is less pronounced. 

(ii) the intensity line ratios will have significant differ- 
ences for both molecules. 

Therefore, the densities derived from astronomical ob- 
servations of HNC and HCN could be different. Moreover, 
in the optically thin case the lines from these species could 



trace different areas of the cloud. Our calculations indicate 
that the j = 1 — line of HNC could show weak maser ef- 
fects for moderate column densities (low line opacities) and 
densities between 10 5 and 10 6 cm -3 for temperatures above 
30 K (i.e., similar to the behaviour of the HCN j = 1 — 
line). However, no population inversion, such as the one ten- 
tatively reported by Aalto et al. (2009) in the direction of 
Arp220, has been found for the j = 3 — 2 line. 

In order to quantify the effect of the new rate coef- 
ficients on the excitation of the rotational levels of HNC, 
we have compared the line intensities obtained with the 
HNC-He rates and those obtained when the HCN-He rates 
are adopted for HNC. Fig. 8 shows the intensity ratio 
R = T B (j 1)[HNC - He]/T fl (j -»■ j - 1)[HCN - He] 

for the j = 1 — 0, 2 — 1, and 3 — 2 transitions. Three differ- 
ent column densities have been considered: 10 12 , 10 13 , and 
10 14 cm -2 (optically thin, r ~1, and optically thick respec- 
tively). The volume density has been selected for the case 
of low, medium and high collisional excitation (n(H2) = 10 4 , 
n(H 2 ) = 10 5 , and n(H 2 ) = 10 6 cm -3 respectively). We see that 
the intensity ratios predicted by both sets of collisional rates 
approach unity for high temperature for all column densities. 
However, for low and moderate kinetic temperature there 
are significant differences in the derived brightness temper- 
atures. The j = 1 — transition is predicted up to 2.5 times 
stronger using the HNC-He rate coefficients than that ob- 
tained from the HCN-He ones. The effect is less prominent 
for the j = 2 — 1 transition and never exceeds 10% for the 
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j — 3 — 2 one. This behaviour has strong implications for 
the interpretation of the observations of HNC in cold dark 
clouds as the density required to obtain the same bright- 
ness temperature is lower for the HNC-He rates than for the 
HCN-He ones. Moreover, using the latter rates it will be diffi- 
cult to fit the line intensity ratios R 12 = T B (\-0)/T B (2-\), 
Ris =T s (l-0)/T s (3-2) and R 23 = T B {2 - 1)/T S (3 - 2). 
R12 has a value of ^ 7 for all column densities, n(H2)^10 4 
cm" 2 , and Tx = 10 K, using the HNC-He rates. Its value is 
~ 5.5 using the HCN-He ones. R23 reachs a value of 17 with 
the HNC-He rates and of 14 with the HCN-He ones for the 
same range of physical conditions. For n(H2) larger than a 
few 10 6 cm" 3 the value of R jr [HNC-He] /R jf [HCN-He] ap- 
proaches unity. For Tk 30 K the R23 ratio does not show 
any significant difference (below 20%) between both set of 
collisional rates. However, the R12 and R13 ratios still show 
important discrepancies for n(¥L 2 ) < 10 5 cm -3 . 

As an example to quantify how important are these 
errors one could consider a typical dark cloud with 
n(H 2 ) = 10 4 cm" 3 , T K = W K and iV(HNC) = 10 13 cm" 2 . To 
fit the predicted j = 1 — line intensity from the HNC- 
He rates using the HCN-He ones we will need to increase 
the column density of HNC by a factor 2.4. Consequently, 
it seems obvious that the HNC-He collisional rates have to 
be used over the HCN-He ones at least for low and moder- 
ate temperatures. Moreover, for low kinetic temperature the 
abundances of HNC reported in the litterature have to be 
revised as the new rates indicated that the HNC abundance 
has probably been overestimated by a factor 2-3 depending 
on the physical conditions of the cloud. It is particularly in- 
teresting to note that the j = 1 — line of HNC, being more 
easily excited than that of HCN, could be formed in zones 
of the cloud different from those of HCN. 

The large different collisional excitation between the 
,7 = 1 — line of HNC and HCN has important consequences 
for the chemical modeling of dark clouds. Irvine & Schloerb 
(1984) pointed out the need for a HNC/HCN abundance ra- 
tio > 1 in these objects. Many other authors have arrived 
to the same conclusion. However, from a detailed model- 
ing Herbst, Terzieva & Talbi (2000) concluded that, from a 
chemical point of view, it was not possible to reproduce the 
derived abundances in dark clouds and that the HNC/HCN 
abundance ratio could not be very different from unity at 
Tk = 10K. The present results provide an explanation to this 
observational problem and allow us to conclude that in cold 
dark clouds the HNC/HCN abundance ratio is ~l. 

For T K > 100 K, it appears from Fig. 8 that the HCN- 
He collisional rates will produce reasonable results in inter- 
preting the observed HNC line intensities and line intensity 
ratios. Hence, the relative abundances between both species 
derived in the literature toward warm molecular clouds are 
not significantly affected by the use of these rates. 



6 SUMMARY AND CONCLUSION 

We have used quantum scattering calculations to investigate 
rotational energy transfer in collisions of HCN and HNC 
molecules with He atoms. The calculations are based on 
new, highly accurate 2D potential energy surfaces. Rate co- 
efficients for transitions involving the lowest 8 levels of these 
molecules were determined for temperatures ranging from 5 



to 100 K. Strong propensity rule for even Aj were found in 
the case of HCN-He system whereas a propensity rules for 
odd Aj were found in the case of HNC-He system. 

The impact of these new rate coefficients on astrophys- 
ical modeling has been evaluated. The new HCN-He colli- 
sional rates have little impact on the interpretation of ob- 
servations for warm and hot molecular clouds as the use of 
the previous rates (Green &: Thaddeus 1974) produce line 
intensities that differ from those obtained from our rates 
by an amount similar to the calibration accuracy of radio 
observations in the millimeter and submillimeter domains. 
However, for low kinetic temperature the previous collisional 
rates could not reproduce the observed line intensity ratios 
j = 1 - 0/j = 2 - 1 and j = 1 - 0/j = 3-2. For Tk<100 K 
we recommend to use our rates. The situation is different for 
HNC as no previous collisional rates are available and the 
interpretation of astronomical observations has been based 
on the use of the Green & Thaddeus (1974) rates for HCN. 
We have found significant variations in the predicted line 
intensities and line intensity ratios using the HNC-He rates 
instead of the HCN-He ones (our calculations for HCN). 
In particular, in cold dark clouds the interpretation of the 
HNC line intensities will dramatically suffer if the HCN-He 
rates are used rather than the correct HNC-He ones. The 
main collider in molecular clouds is however molecular hy- 
drogen. Hence, these conclusions could be limited to the case 
of HCN/HNC-para-H2(j = 0) for which the collisional rates 
are expected to be very similar to those of HCN/HNC-He 
(see below). 

Our results for HNC indicate that the HNC/HCN abun- 
dance ratio derived from observations toward dark clouds 
has to be revised. The best fit to the observations available 
in the literature is obtained for HNC/HCN ~ 1. We have 
neglected in our calculations the hyperfine structure of HCN 
and HNC. Although difficult to spectroscopically resolve for 
the latter, for the former species the line intensities of the 
different hyperfine components are known to show impor- 
tant intensity anomalies (Kwan & Scoville 1975; Walmsley 
et al. 1982; Cernicharo et al. 1984; Gonzalez- Alfonso & Cer- 
nicharo 1993), which can be explained as the result of their 
different opacities across the cloud (Cernicharo et al. 1984; 
Cernicharo & Guelin 1987; Gonzalez- Alfonso & Cernicharo 
1993). Line overlaps also play an important role in the ex- 
citation of HCN in molecular clouds (Gonzalez- Alfonso & 
Cernicharo 1993). They are responsible of the strong masers 
found in other species (see, e.g., Gonzalez- Alfonso et al. 
(1996); Gonzalez-Alfonso & Cernicharo (1997)). Hence, a 
correct treatment of the radiative transfer of HCN (and in 
some cases also of HNC) in molecular clouds will require a 
detailed knowledge of the collisional rates between hyperfine 
components. These calculations are under development and 
will be published elsewhere. 

Finally, as already said, the great abundance of H2 in 
the interstellar medium makes this molecule the primary col- 
lision partner for any other species. It is generally assumed 
(Lique et al. 2008) that rate coefficients with He can provide 
a good estimate of rate coefficients for collision with para- 
H2O = 0). This approximation postulates that collisional 
cross sections with He and para-H2(j=0) are equal, so that 
the rate coefficients differ only by a reduced mass factor of 
« 1.4 arising from the thermal averaging (Eq. 3). Recent 
results on rotational excitation of CO (Wernli et al. 2006), 
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SO (Lique et al. 2007) and SiS (Lique & Klos 2008; Klos & 
Lique 2008) have pointed out that rate coefficients for colli- 
sions with para-H2 (j = 0) can be up to a factor of 3 larger or 
smaller than those for collisions with He, depending on the 
selected transition, but that the He rate coefficients scaled 
by a factor f.4 provide the correct order of magnitude of 
the H2C7 = 0) rate coefficients 4 . Therefore, the present re- 
sults should provide a reasonable first estimate of collisional 
rate coefficients for collisions of HCN and HNC with para- 
H2O' = 0). On the other hand, it may be unadvisable to 
use the present He rate coefficients as an estimate of the 
ortho-H2 rate coefficients, since the He and ortho-H2 rate 
coefficients usually differ significantly. Specific calculations 
with ortho-H2 must be performed. Thus, it is crucial to ex- 
tend the calculations, both of the PES and of the inelastic 
cross sections, to the HCN-H2 and HNC-H2 systems. 

In any case, the full set of rate coefficients presented 
here will enhance our ability to understand and interpret 
future HCN and HNC observations. Collisonal rates when 
the hyperfine structure is considered can be obtained from 
the present results and the standard assumption that the 
rates will be proportional to the degeneracy of the final hy- 
perfine state, or employing the more physical IOS scaling 
method (Faure et al. 2007). 
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